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Optical Schemes for Quantum Computation in Quantum Dot Molecules 
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We give three methods for entanghng quantum states in quantum dots. We do this by showing how 
to tailor the resonant energy (Forster-Dexter) transfer mechanisms and the biexciton binding energy 
in a quantum dot molecule. We calculate the magnitude of these two electrostatic interactions as 
a function of dot size, interdot separation, material composition, confinement potential and applied 
electric field by using an envelope function approximation in a two-cuboid dot molecule. In the 
first implementation, we show that it is desirable to suppress the Forster coupling and to create 
entanglement by using the biexciton energy alone. We show how to perform universal quantum logic 
in a second implementation which uses the biexciton energy together with appropriately tuned laser 
pulses: by selecting appropriate materials parameters high fidelity logic can be achieved. The third 
implementation proposes generating quantum entanglement by switching the Forster interaction 
itself. We show that the energy transfer can be fast enough in certain dot structures that switching 
can occur on a timescale which is much less than the typical decoherence times. 



PACS numbers: 03.67.Lx, 03.67-a, 78.67.Hc, 73.20.Mf 



I. INTRODUCTION 

Quantum dotgi^ are quantum heterostructures which 
are composed of nanoscale regions of one type of ma- 
terial which is embedded in a second type. In a semi- 
conductor quantum dot (QD), materials with differing 
bandgaps are used; this leads to the possibility of elec- 
tronic confinement within the dot region. Moreover, the 
confined electronic states can be accurately controlled by 
varying the dot size, shape or composition, and the num- 
ber of confined electrons; all of these may be altered by 
using different growth conditions and hence specifically 
tailored "artificial atoms" or "superatoms" can be pro- 
ducedi^ Some prominent atom-like properties of QDs 
include an electronic shell structure,"^ Rabi oscillations)^ 
photon antibunching,^'^ controlled quantum light emis- 
sion,^'* and quantum entanglement. ^'^^ One of the most 
intriguing possible applications of quantum dots is that 
they may be used to build quantum computers^ii^ A 
practical realization of a quantum computer would be 
very significant, since there exist theoretical quantum al- 
gorithms which would make some classically hard com- 
putational problems tractableii^ Such quantum devices 
could also accurately simulate any physical system (and 
the evolution of its local interactions) by invoking the 
same amount of energy and Hilbert space requirements 
as the system itselfii^*i^*i& 

The basic unit of a quantum computer is a two-level 
quantum system, the so-called qubit. Of the utmost im- 
portance is the identification of a physical system where 
a coherent qubit evolution can be performed, thus al- 
lowing a precise execution of the elementary quantum 
gates required for universal quantum computation. ^^ii^ 
Many different types of hardware for embodying qubits 
have been proposed (for a collection of papers detailing 
some of these see Ref. UTI) and some of them have al- 



ready been implemented for performing elementary quan- 
tum gate operations. These include ion traps fi^ii^iSSii^MS, 
quantum electrodynamics cavitiesr2i2ii2Si2S42i nuclear 
magnetic resonancefSSiSSi^Si^, dopants in semiconduc- 
^Qj.g^32j33j34 optical lattices and Bose-Einstein condcn- 
gg^^gg^35j36j37 Josephson junctions^SSi^SiiMi and quantum 
dotsii°i4^i^3,44,45,46,47,48,49^50,5i ^^vis article we concen- 
trate on a quantum dot implementation. Previous pro- 
posala9d°.42i43,44,46,47,5i,52,53 -^^^x^^q the usc of a sin- 
gle electron^SiSiiS or nuclear— spin located on each of 
an array of interacting dots, the presence or absence 
of an electron charge state^^ or the use of excitonic 
states— i^^i^^i^ We show how an energy selective ap- 
proach to manipulating the excitonic states of coupled 
QDs, together with control over the energy transfer and 
biexciton binding energy, can be used to perform quan- 
tum computation (QC) and to produce controlled exci- 
ton quantum entanglement. In so doing, we investigate 
the Forster-Dexter resonant energy transfer, a mecha- 
nism first studied in the context of the sensitized lu- 
minescence of solids)^^*^ in which an excited sensitizer 
atom can transfer its excitation to a neighbouring ac- 
ceptor atom, via an intermediate virtual photon. This 
mechanism is also responsible for photosynthetic energy 
processes in antenna complexes, biosystems (BSs) that 
harvest sunlight. More recently, interest has focussed 
on energy transfer in quantum dot nanostructures^^ and 
within molecular systems (MSs) — In this article we show 
how to exploit such energy transfer mechanisms with a 
view to processing quantum information. This article is 
a more detailed account of the work which appears in 
Ref. EE 
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II. BUILDING QUANTUM LOGIC GATES 

We consider the Hamiltonian of two interacting quan- 
tum dots. We assume that the dots are sufficiently far 
apart that tunneUing processes between them may be 
neglected but that there is a strong exciton-exciton cou- 
pling. Our two-level system is represented in each dot by 
a single low lying exciton state |1} and the ground state 
0). Then the interaction Hamiltonian can be written in 
the computational basis ({|00}, |01), |10}, |11}}, with the 
first digit referring to dot I and the second to dot II) as 
follows: 
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The diagonal interaction Vxx is the direct Coulomb bind- 
ing energy between two excitons, one located on each 
dot, and Vp denotes the Coulomb exchange (Forster) in- 
teraction which is off-diagonal and therefore induces the 
transfer of an exciton from one QD to the other. These 
are the only Coulomb interaction terms which act be- 
tween the qubits and will be calculated and discussed in 
detail in Section ioq denotes the ground state energy, 
uj\ (W2) refers to the energy required to create an exciton 
on dot I (II) in the absence of interactions, and includes 
intra-dot coupling contributions (direct Coulomb bind- 
ing energy and spin splitting) which we shall discuss in 
Section IIVI We also define Ao = wi — W2 to be the dif- 
ference between the exciton creation energy for dot I and 
that for dot II in the absence of interactions between the 
dots. Thus, if Ho = Hi + Hu denotes the free particle 
Hamiltonian, then -ffo(l7i) I72)) = (7iWi-t-72W2) I71) I72), 
71,72 = 0,1 {h — 1 throughout this article), then 
H — Hf) + Vi^2 (where Vi^2 accounts for the qubit-qubit 
interactions, Vxx and Vp) is the system's overall Hamil- 
tonian. In the case of an n-qubit register with nearest 
neighbour interactions, the Hamiltonian takes the form 
H = if^"^ + where 7?^"^ = EtlH^ is 

the free particle Hamiltonian, and are the inter- 

action terms. A related Hamiltonian was investigated 
in Ref. [s^l, but there the off-diagonal interaction terms 
(Vf) were neglected. 

The eigenenergies and eigenstates of the interacting 
qubit system are 



1*00) = |00); 



£■00 — 

Eqi = Wo H 

Ew = wo + wi - 4^(1 - A), 1*10) = -ci |01> 

£^11 = Wo + wi + W2 + Vxx, 1*11) = 111) , 



Wl-^(1+A), 1*01) =Cl|10)+C2|01) 



C2IIO); 
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where A = I + 4{Vf / A^y , a = ^(A - 1)/2A 
(w Vf/Ao for Vf/Ao < 1) and C2 = ^/{A + 1)/2A . The 
eigenenergies in the absence and presence of interdot 
interactions are displayed in Figs. H^a), (b) and (c). 
Fig. n^a) shows the energy levels when the interactions 



are off; Fig. Q^b) shows these when the interactions are 
on, but where Vf <^ Aq; Fig.^c) shows Eiq and £^01 as 
a function of the ratio Vf/Aq. Fig. ^d) shows ci and C2 
as a function of Vf/Aq. These figures demonstrate that 
Vp causes a mixing of the states |01) and |10) such that 
the eigenstates of the interacting system are not the 
same as the computational basis. As we show below, this 
Vp coupling can be used for generating highly entangled 
states. 

Single qubit operations can be achieved by inducing 
Rabi oscillations in the excitonic system (e.g., see Refs.|j 
and l5^ . If we take a Bloch sphere representation of 
a qubit, where the state |0) is represented by a unit 
vector from the origin to the north pole of the Bloch 
sphere and the state |1) by a unit vector to the south 
pole, then the qubit state \tp) = exp(iA)(cos(0/2) |0} -I- 
exp(z(/?) sin(6'/2) |1)), where A is a global phase, is defined 
by the unit vector (cos(/3sin6',sin(^sin0,cos0). Depend- 
ing of the values of 9 and this vector may point to 
any point on the surface of the sphere, and for univer- 
sal quantum computation it must be possible to move 
the vector between any two of these points; this defines 
an arbitrary single qubit rotation. In our QD system, 
this control can be achieved by using laser pulses to in- 
duce two distinct Rabi oscillations (see Fig. The en- 
ergy and length of such pulses must take into account 
structural factors like the dot confinement energies and 
transition dipole momentsi^2iSS It is also essential that 
the exciton states have long enough decoherence times 
that control over the phase tp is possiblei^i*^ For exam- 
ple, self-assembled semiconductor (e.g. InCaAs/GaAs) 
quantum dots could be advantageous for qubit manipu- 
lations since they exhibit large dipole moments and long 
dephasing times. ^'^^ We shall return to the role of the 
QD material composition parameters for QC below. 

The Vxx and Vp interactions lead to three possible 
ways of achieving quantum entanglement. First, if the 
ratio Vf/Aq ^ 1, the eigenstates of the system are ap- 
proximately |00) , -l=(|io)-|01)), -i=(|10)+j01)) and |11). 

We now further assume that the ratio Vp / Ao can be con- 
trolled, by means of applying an electric field to either 
change V^f directly, or to increase Ao by means of the 
Stark shift (we shall discuss both of these effects in detail 
in Section fV A|) . Then, we initially prepare the system in 
a state where Vf/Aq <^ 1 and we selectively excite QD 
I and create |10). Now, when the Forster interaction is 
turned on, the system will naturally evolve sequentially 
into the following states: |10) ^ 7f(|10) + « |01)) ^ 
|01) ^ 7f(|10) - «|01)) ^ |10) (see Fig.|g. This evo- 
lution could be then stopped when the system is in a 
maximally entangled state by applying an electric field to 
suppress the Forster coupling once more, an effect which 
we shall again explore in detail in Section IV Al 

Second, if the system does not have a strong Forster 
coupling, i.e. Vf/Aq 1 all the time, the computational 
basis states are essentially the eigenstates of the system. 
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FIG. 1: A schematic diagram of the properties of the model Hamiltonian, Eq. Q (a) energy levels in the absence of qubit 
(interdot) interactions; (b) energy levels in the presence of qubit interactions for two dots (I and II) of different excitation 
frequencies. ei2 = uj2 + Vxx — S, e2i — loi + Vxx + S, and S = Vj?/Ao, where Vf and Vxx represent the strength of the Forster 
and direct Coulomb binding interactions respectively. In this case Vf <^ Aq. (c) Eigenenergies -Boi, and -Bio corresponding to 
the qubit eigenstates l^'oi), and l^'io) as a function of the ratio Vf/Aq for wi/Ao = 20. For comparison, the dashed lines show 
the energies when Vf/Aq — 0. (d) The eigenstate coefficients Ci as a function of the strength Vf/Aq. These coefiicents show the 
departure from the basis states |01), and |10) followed by the eigenstates j^oi) = ci |10) +C2 jOl) and j^io) = — ci |01) +C2 jlO). 
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FIG. 2: Bloch sphere representation of a single qubit. In 
order to perform an arbitrary single qubit rotation in the ex- 
citon system, it is necessary to have the ability to induce two 
different Rabi oscillations. Example trajectories for two such 
oscillations are shown by dashed lines. 



Then the Vxx couphng term implies that the resonant 
frequency for transitions between the basis states |0) and 

1) of one qubit depends on the state of the neighbouring 
qubit. This means that it is possible to construct a CNOT 
gate in this system. Such a gate flips the target qubit k 
if the control qubit j is in the state |1) and acts trivially 
otherwise: CNOTjk{\m)-\n)i^) i-^ \m)-\m®n)j^, where 
Tn,n ^ {0, 1}, and © denotes addition modulo 2 or XOR 
operation. By referring to Fig. Q^b), we can see that 
the logic operation CNOTi2(|1};^ 10)2) ^-^ |1);^|1)2 can be 
achieved by illuminating the qubit system in the state 

10) with a TT-pulse of energy £12 — UJ2 + Vxx — 5 {a pulse 
which we label tt^jj). Conversely, if the role of the control 
qubit is to be performed by the second qubit, the gate 
operation CNOT2i(|0)j 11)2) 1— > |1)2 can be reahzed 
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FIG. 3: Illustration of how to create an entangled state by 
manipulating the off-diagonal Forster coupling between two 
nanostructures. When the interaction is on and much greater 
than Ao , the eigenstates of the system correspond to the large 
dots on the equator of the Bloch sphere. Thus after selective 
excitation of state 1 10) the system will naturally evolve to the 
equator. The interaction may then be suppressed by means of 
an applied field, and the state remains maximally entangled 
as the eigenstates become the computational basis states. 



via the application of a 7r-pulse of energy £21 = uji+ Vxx + 
S (a TTe^j pulse) to the system state |01). Crucially, the 
energy difference between these two 7r-pulses, and hence 
the energy selectivity of the logic gate, is determined by 
Ao, Vxx and Vp- From Fig.^b) we can also see how to 
use the CNOT gate to create maximally entangled states. 
For example, if we start in the ground state and first 
apply a n/2 or 37r/2 pulse at energy wi, we create the 
states --i=(|00) ± |10)); if we now apply a n^^^ pulse, we 

generate the maximally entangled states -i=(|00) ± |11)). 
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Third, the interaction with the laser field does not 
necessarily have to be the entangling mechanism in the 
case where Vf/Aq <C 1. If each of two neighbour- 
ing single qubits are prepared in the superposition state 
-^(|0) + |1)) (i.e. a Hadamard transform is applied to 

both of them) , making the state i ( 1 00) + 1 01 ) + 1 1 0) + 1 1 1 ) ) , 
this will then naturally evolve into entangled states under 
the action of Vxx alone. 

Although we have discussed our Hamiltonian specifi- 
cally in the context of a dot molecule, these ideas are 
valid for any nanostructure with corresponding diagonal 
and off-diagonal interactions where an energy selective 
excitation is possible. However, we focus on a quantum 
dot implementation, and in the next sections we shall an- 
alyze in detail the inter- and intra-dot interaction terms 
and show how they can be tailored as a function of the 
interdot distance, dot sizes, and material composition. 



III. THE MODEL SYSTEM 

Different dot geometries (e.g. spherical, pyramidal or 
cuboidal shaped dots) can be used to implement the 
logic gates and quantum entanglement schemes discussed 
above, and in this article we choose dots of a square- 
based cuboidal shape (for another possible geometry see 
Ref. 62). We assume that the potential energy V of 
both electrons and holes increases abruptly at the cuboid 
boundaries where the semiconductor bandgap changes, 
and that V = inside the cuboids (see Fig. QJ- The 
confinement potential is determined by the band offsets 
for the electrons and the holes. This type of square well 
potential has the advantage of describing both a well de- 
fined dot size in all three dimensions and of bound and 
unbound solutions in each direction (this is in contrast 
with, for example, the parabolic potential considered in 
Ref.lHl). 

Our model captures the essential properties of quan- 
tum dots which are grown by the Stranski-Krastanow 
(SK) methodiS^ Such structures show a degree of self- 
organization^^ which is ideally suited for the manufac- 
ture of prototype quantum devices, and in the realization 
of the elementary quantum logic gates discussed in this 
paper. The SK dot growth proceeds through the evapo- 
ration of a layer of dot material on to a previously grown 
substrate: dots form spontaneously due to the competing 
energy contributions of dot surface area, dot volume and 
strain during the growth process. After the formation of 
dots and subsequent overgrowth of the substrate mate- 
rial, a second layer of dots may be grown; these nucleate 
preferentially above the dots in the first layer due to the 
uneven strain field at the surface. An example of two 
dots grown in this way is shown in Fig. [S] ~ and it can 
be seen there that two vertically stacked dots have been 
grown with a controlled spacing; dots with such char- 
acteristics may be well suited for performing the logic 
operations described above. As can be seen in Fig.|Sl the 
dots tend to have smaller dimensions in the growth direc- 
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FIG. 4: Schematic diagram of the cuboidal dot model. The 
cuboids have base sides of length 2a and 2b and heights of 
hi and h2 respectively. Their centres are separated by a dis- 
tance R. The potential inside the cuboids is set to zero, and 
that outside them is determined by the band offsets of the 
conduction and valence bands within the heterostructure. 




FIG. 5: Transmission electron micrograph of two layers of 
quantum dots grown by the Stranski-Krastanow method. The 
dots are made of InAs and the encapsulating material is GaAs; 
note that the dots in the second layer nucleate preferentially 
above the dots in the first layer. (Figure courtesy of C. Lang, 
C. Marsh and D. Cockayne (Oxford Materials); sample cour- 
tesy of M. Hopkinson and P. Houston (Sheffield University).) 



tion than in the perpendicular directions, and the upper 
dot of the pair tends to be of a slightly larger size, thus 
allowing for an appropriate identification of the excita- 
tion frequencies ui and UJ2 required in our model. Qubit 
scalability is available via the SK growth procedure since 
several layers of dots have been shown to grow in stacksiS^ 
Our computational Hilbert space requires up to one ex- 
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citon per dot, and we must therefore calculate both the 
single particle energies, which are determined by the po- 
tential profile of the dots, and the two- and four-particle 
interactions which are determined by the strength of the 
Coulomb interaction between them. We look at these 
different quantities separately in this section. 



A. Single Particle States 

There are a variety of methods for finding the solu- 
tions of the Schrodinger equation for electrons or holes 
in a quantum dot. These include a full pseudopotential 
calculation(S£iSS*SS finite element analysisjBi plane wave 
expansion'''^ and the use of finite differences We em- 
ploy a similar strategy to that in Refs. and Fl. where 
the Schrodinger equation is expanded in a set of ana- 
lytical basis functions which are the exact solutions of a 
potential which is close to the one under investigation. 
This method has the advantage that the state solutions 
can be stored in a vector in Hilbcrt space rather than as 
a wavefunction amplitude at each of a very large number 
of different spatial points; time evolution of the quan- 
tum states is also easier to simulate when the state is 
represented by a vector and we shall extend our work to 
this area elsewhere. Furthermore, a vector representation 
can allow more physical insight since the basis functions 
themselves have a known physical interpretation. 

Our first step is to express the wavefunctions for single 
particles in the envelope function approximation asi2^ 



a=-Sh I 



V'p(r) = (/)p(r)C/p(r) 



(3) 



where 0p(r) is an envelope function describing the chang- 
ing wavefunction amplitude of confined states for particle 
type p over the dot region, and Up{r) is the Bloch func- 
tion which has the periodicity of the atomic lattice. In 
the effective mass approximation, the envelope functions 
are solutions of the following single particle Schrodinger 
equation: 



/ 1 
[m;{r) 



V + Vp{r) 



^;{r) ^ E;^;ir) (4) 



where Vp is the confinement potential, which is displayed 
in Fig. 01 and m* is the effective mass of the particle 
p. These solutions may be obtained by expanding the 
Hamiltonian in a set of envelope basis functions of the 
form S(r) — ^x{x)iy{y)£,z{z), where the are the so- 
lutions of a one dimensional square well potential with 
the appropriate effective massesi^SiZi Both bound and 
unbound states must be used in the expansion in order 
to obtain convergent solutions: the forms of these are 
discussed in Appendix A. 

There are two important things to mention about the 
direct expansion technique. First, the basis functions we 
have described above do not in general form an orthog- 
onal set if the bound state solutions have not decayed 
to zero at the artificial infinite barrier which is used to 
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FIG. 6; The single particle ground state energy as a function 
of dots size for a quantum cube (lower graph) and a quan- 
tum cuboid (upper graph). We use the envelope function 
and effective mass approximations. The filled symbols cor- 
respond to nip = 0.6mo which is typical of heavy holes, and 
the open symbols correspond to rrip = 0.06mo which is typ- 
ical of electrons. These values will be used for electrons and 
holes throughout this paper. The different symbols represent 
difi'erent values for the confinement potential (see legend). 



generate the unbound states. In practice this is rarely a 
problem, but anyway is circumvented by using a modi- 
fied basis set which is orthogonal and whose components 
are linear combinations of the original basis functions. 
This modified basis set spans the same Hilbert space as 
the original set. The method we employ to find this set is 
canonical orthogonalization^^ which relies on direct diag- 
onalization of the matrix whose elements are the overlap 
integrals of the basis functions. The Hamiltonian may 
then be expressed in this new basis as an Hermitian ma- 
trix, and solutions arc found again by direct diagonaliza- 
tion (we use the NAG diagonalization algorithm in our 
simulations) . The second point is that the basis set must 
necessarily be truncated; hence the eigenenergies of the 
solutions we obtain are really upper limits on the true 
eigenenergies of the coupled dot system (we employ the 
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Rayleigh-Ritz variational methoc'J^). In practice we can 
increase the number of basis functions until a sufficiently 
accurate solution is obtained. We shall mention any im- 
portant points relating to solution convergence and ap- 
proximations at the appropriate places later in the paper. 

Results of a simple single particle calculation are dis- 
played in Fig. where the ground state energy of a two 
particles in QD I (of masses 0.6mo and 0.06mo, where toq 
is the free electron mass) is shown as a function of the dot 
size for two dot geometries. The two geometries corre- 
spond to a cubic shape (a = hi/ 2) and a flatter cuboidal 
shape (a = 5/ii), which is more typical of SK dots (see 
Fig. |3J). As would be expected, the ground state energy 
decreases for a larger dot size and is smaller for heav- 
ier particles. The cube shaped dots have ground states 
with smaller energies than the corresponding cuboidal 
dot states since the cuboids have one smaller dimension 
which increses the kinetic energy of the wavefunction. All 
of the curves have a kink, and for dot sizes below this the 
ground energy saturates to the value of the confinement 
potential. This is a consequence of including unbound 
states in the calculation: once the dots are small enough 
that the confined state energy has a larger energy than 
the confinement potential, the ground state becomes un- 
bound and there is hardly any dependence on dot size as 
it is further reduced. When the ground state is a bound 
state, the ground energy is always larger in the case of 
the cuboidal shape, for a given value of a: this is a con- 
sequence of the smaller height dimension of the cuboid 
which increases the kinetic energy of the particle. Bound 
ground states are always very closely approximated by 
the basis function corresponding to the ground state of 
the one dimensional well in all three dimensions: in the 
case of cubic shaped dots, the amplitude of this state is 
always greater than 0.999; for the cuboidal dot it is al- 
ways greater than 0.99. We shall use the approximation 
that the ground state is exactly this basis function later 
in the paper. 



B. Coulomb Interactions: General Methodology 

In this section we present our general methodology for 
calculating the Coulomb interaction matrix elements be- 
tween electrons and holes in quantum dots. This will be 
important later for calculations of the intra- and interdot 
matrix elements in our model system. 

First, consider an initial wavefunction of an iV-electron 
system which represents a single exciton state of a quan- 
tum dot:*^^ 



V'l(ri,cri),V'2(r2,cr2), 

V't(rt,crt), ...,'(/'Ar(rN,crN) , (5) 



where the A indicates overall antisymmetry (i.e., the 
wavefunction takes a Slater determinant form), the Ui 
represent the spin state of each electron, and the ipi are 



single particle wavefunctions; we have labelled the state 
s with a prime symbol to indicate that it lies in the con- 
duction band, whereas all of the other states are in the 
valence band. 

Next, we assume a final state which is a different single 
exciton: 

^-F = A i>i(Ti,ai), 11:2(^2, 02),-; 

■0s(rs,crs), ...,i/'J(rt,crt), ...,?/'Jv(rN,crN) . (6) 

The Coulomb matrix element between these two states 
is given by: 



ij-i<j 



47reoer(rij)|rij 



1*1), (7) 



where eq is the permittivity of vacuum, er(r) is the rela- 
tive permittivity of the medium (and therefore describes 
polarization screening), and rjj = r; — rj. The only non- 
zero terms in the above expansion are those involving 
both Fg and rt, since the ground and excited states of 
each single particle are orthogonal to one another. Hence, 
we obtain: 



47reoer(rst)|rst| 



1*1) 



(8) 



Owing to the antisymmetric nature of the wavefunctions, 
this matrix element has contributions from a direct term 
and an exchange term. Both terms take the form of an 
integral over Ys and rt. In spite of the fact that it arises 
from the exchanged form of the wavefunction, the direct 
term is conventionally written as: 



Mii = C//C(r.)^;(rs) 



1 



er(rst) 



I?/'* (rt)-0t(rt)rfrsdrt, 
(9) 



and the exchange term is 



±c hiij^M^s) 



er(rst)|rst 



(10) 



where we have introduced the constant C — — . The 

47reo 

sign of this exchange term is determined by the spin of 
the two particles: spin triplet [S = 1) states have positive 
Mj^ elements, whereas spin singlet (S = 0) states have 
negative values. We have removed the spin vari- 
ables from the single particle wavefunctions since these 
are not important for spatial integrals, once the sign for 
Mj^ has been determined. It is natural now to switch to 
a hole description of the many body wavefunctions intro- 
duced above: the matrix elements obtained by the above 
procedure would be identical if we labelled the states s 
and t as holes and simply expressed each wavefunction 
as a product of the promoted electron states and the left 
behind hole states. This allows the other parts of the 
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wavefunctions to be left out when calculating matrix el- 
ements involving these specific electronic states and so 
makes calculations easier; it also explains the convention 
on labelling direct and exchange terms. We shall hence- 
forth use such a description. 

We proceed further by noting that the two integrals, 
Eq. 1^1 and Eq. ^| are both of the form 

I = JJ Ps{rs)f{rs - rt)ptirt)drsdrt, (11) 

and so we may use Fourier transforms to reduce the di- 
mensionality of the integrand. Employing the convolu- 
tion theorem and Parseval's relation leads to 

I = J RsiK)R,(K)F(K)dK, (12) 

where the Fourier transform of p is denoted by R and 
that of / is denoted by F. We now make the calculation 
more specific by first assuming that is independent 
of Tst (see Ref. jG^j for a detailed discussion of the form 
of er(rst)), and note that the Fourier transform of the 
Coulomb operator is given by 

F{K) = 1^ . (13) 

The p functions are a product of two wavefunctions of the 
form of Eq.O and their Fourier transform is simplified by 
invoking the different lengthscales of the envelope func- 
tion and the Bloch function. We may write, for the Mjp 
integral (a completely analogous method can be carried 
out for the Mjp integrals): 

i?.(K) = y"v:(rs)V'Urs)exp(iK-r,)drs 
= Vceii ^ rsA,^ exp(jK • TO X 

Ti 

/ C/;(r,)C/^(r,)exp(zK.r,)dr, , (14) 

Jcell 

where we have assumed that each envelope function takes 
a constant value, 0j, over each unit cell i of the lattice, 
that the translational lattice vector of cell i is and 
that the volume of the unit cell is Vceii- We may express 
the K wavevector as 

K = k + G, (15) 

where k is a vector within the first Brillouin zone and G 
is a reciprocal lattice vector, and we may further convert 
the sum to an integral to obtain 

Rs{k,G) = I (j)l{r)(j)',{r) exp{ik ■ r)dr X (16) 

•/ space 

f C/;(r,)[/^(r,)exp(z(k-f G) •r,)dr, . 

Jcell 

Thus Rs is a product of Fourier transforms, one for the 
envelope functions (which is independent of G, so only 



needs to be calculated within the first Brillouin zone) 
and the other for the Bloch functions. It is also obvious 
that an analogous expression exists for Rt . The envelope 
function parts of the Fourier transform are analytical: 
the wavefunction takes either a sinusoidal or exponential 
form depending on whether it exists within or outside 
the quantum dot. We do not write out these expressions 
explicitly here however, since they are somewhat lengthy 
and tedious. The Bloch function part is also analytical 
for a suitable choice of wavefunction, and to simplify the 
calculations here we use a Kronig-Penney model where 
the atomic wavefunctions are assumed to take the form 
of the solutions of an infinite square well potential of well 
width 2x. Specifically, we assume that the hole states we 
consider (at the top of the valence band) take the wave- 
function solution of this potential which has Pz symmetry 
and the electron states we consider take the solution with 
s symmetry. We expect this to be sufficient approxima- 
tion for elucidating the general properties of the system, 
though a more refined calculation would be required to 
obtain more accurate estimates of the various quantities 
we calculate. 

By inserting Eq. 1161 (and the analogous expression for 
Rt) and Eq. El into Eq. E|wc obtain an expression for 
Mj^. The expression has an integrand which is analyt- 
ical but the integration over three-dimensional K space 
must be carried out numerically. This is done by employ- 
ing a NAG library routine for multi-dimensional adap- 
tive quadrature. If the wavefunction labels are swapped 
around, an analogous method for calculating Mjp can 
be carried out (see Eq. EJ. We shall show the results of 
such calculations for various different cases in the next 
sections. 



IV. INTRA-DOT COUPLING 

In this section, we describe the predictions of the above 
model when it is applied specifically to the calculation of 
the diagonal matrix element of the two ground state ba- 
sis functions representing an electron and a hole on the 
same dot. In this case the states s and t are identical, 
and hence the expression for Mjp reduces to the direct 
Coulomb interaction between the ground basis state elec- 
tron and the ground basis state hole (we call this Mqq). 
We saw in the previous section that the ground basis state 
is a good approximation to the true ground state of the 
system when only the single particle contributions to the 
Hamiltonian are taken into account. Thus this matrix el- 
ement is a first order correction to the energy due to the 
Coulomb force between the two particles. Furthermore, 
the expression for reduces to the spin splitting be- 
tween singlet and triplet exciton states in this first order 
approximation. 

Let us first consider the direct Coulomb interaction in 
dot I, which has basal side length a, and simplify things 
by calculating for a cubic shape (i.e., we set 2o = hi). 
The results may then be directly carried over to dot II, 
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FIG. 7: The intra-dot direct Coulomb interaction strength in- 
tegrand, plotted in the K^ — Kz plane of reciprocal space, and 
around the K = point. We have calculated the integrand 
for dot I and used the cubic geometry (a = hi/2 = 10 nm), 
and assumed that Ve = Vh = 500 meV. 
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FIG. 8: The intra-dot direct Coulomb interaction strength, 
Mqq as a function of dot size and confinement potential. The 
solid lines are for a cubic geometry (a = fti/2), and the dotted 
lines are for a cuboid {a — 5hi). 



with a ^ b and 2b = h2- As we described earlier, it is nec- 
essary to evaluate the integral of Eq.|^by first transform- 
ing to reciprocal space and then integrating over K space. 
The resultant integrand has peaks at each reciprocal lat- 
tice point (where k = and so the envelope function part 
of the integrand has a maximum). These peaks quickly 
die away over a lengthscale ~ 1/a as would be expected 
for envelope functions representing wavefunctions within 
QDs of side length a. However, it turns out that only 
the central (G = 0) peak is important, since the other 
peaks contribute much less to the total integral (this is 
caused by both the k dependence of the atomic contribu- 
tion to the integrand and by the 1/fc^ dependence of the 
Coulomb interaction part). The central peak is displayed 
in Fig.CI as a function of and Ky [K^ = 0). 

By numerically integrating the central peak for a range 
of dot sizes and confinement potentials, we can obtain 
a plot of the dependence of Mqq on these parameters 
(this is shown in Fig. [S] for both the cubic geometry and 
for a flat cuboid, in which a = 5/ii). As would be ex- 
pected, the interaction decreases as the size of the QD 
increases (and so the electron and hole are not forced to 
be so close together). It is interesting to look at what 
happens at shorter distances when the confinement po- 
tential changes; a larger confinement potential causes a 
larger Coulomb binding energy. This result is expected 
since the wavefunction of both the electron and the hole 
is contracted when the confinement potential is larger - 
and the resultant closer proximity of the two wavefunc- 
tions causes a larger Coulomb interaction. As would be 



expected intuitively, at very large dot sizes the size of 
the Coulomb interaction scales like 1/a. At very small 
dot sizes, the direct Coulomb interaction does not fol- 
low these simple rules: for a weak enough confinement 
there is a peak in the energy and at small values of a 
it decreases. This can be understood by thinking about 
the shapes of the wavefunctions in this region. When the 
well width is small, the curvature of the wavefunction 
is necessarily rather high, and so the kinetic energy is 
large. In order to compensate for this, the wavefunction 
spreads out into the barriers at the cost of some poten- 
tial energy, if the barriers are not too high (this energy 
cost is balanced by the saving in kinetic energy). Thus 
the wavefunction has a larger size than would naively 
be expected from the dot size, and the Coulomb binding 
energy decreases. The peak in the Coulomb potential 
occurs for a larger value of a in the cases of cuboidal 
geometry: this is simply because of the shorter height 
dimension in this case, which means that the wavefunc- 
tion spread effect discussed above remains significant at 
larger values of a. At still larger values of a, the Coulomb 
interaction is larger for the cuboidal geometry than it is 
for the cube: again this is because the cuboid has one 
smaller dimension, which means that the electron and 
hole are forced to be closer together in the cuboidal case. 

We next look at the value of the exchange coupling 
between the two ground state basis functions (i.e. the 
spin singlet-triplet splitting, M^, correct to first order). 
The relevant if -dependent integrand takes a somewhat 
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different form in this case. The central peak (around 
G = 0) is displayed in Fig. O where it can be seen that 
the function has a zero at K = 0; this is expected since 
the electron and hole wavefunctions have opposite par- 
ity. The suppression at if = means that, this time, 
the regions around other reciprocal lattice points have to 
be included in the numerical integration. The resultant 
dependence of Af,^ on dot size and confinement potential 
is displayed in Fig. ^1 

By reference to Fig.llUlwe see that the exchange split- 
ting is several orders of magnitude smaller than the direct 
Coulomb term, though it follows the same trends of in- 
creasing in value with smaller dot sizes and larger confine- 
ment potentials. These effects can be understood as fol- 
lows: the exchange splitting is essentially a consequence 
of Pauli's exclusion principle which states that particles 
in the same quantum state cannot exist together at the 
same spatial point. Thus electrons and holes which are 
in a triplet spin state (and so have indistinguishable spin 
properties) necessarily 'avoid' each other, thus reducing 
the Coulomb attraction between them. This effect is ex- 
pected to be more significant when the wavefunctions in 
the absence of the Coulomb interaction overlap strongly 
- that is when they are localized in one small region of 
space. There is, in general, a greater degree of localiza- 
tion of the wavefunctions when either the dots are smaller 
or when the confinement is stronger, and hence the ex- 
change splitting gets larger when these conditions are sat- 
isfied. However, there is again one region of Fig. [TUl where 
this general rule is not obeyed - that is when confinement 
is relatively weak, but where the dot size is small. Here, 
the exchange energy takes a down turn as the dot size gets 
smaller. This is caused by increased wavef unction barrier 
penetration, which again means the effective wavefunc- 
tion size gets larger rather than smaller as might be ex- 
pected from the simple intuitive picture described above. 
The comparison between the cubic and cuboidal geome- 
try show a similar trend to that discussed above for the 
case of the direct Coulomb interaction, and the reasons 
for this follow the same lines. The turnover occurs for 
a larger a for the cuboidal dot due to the wavefunction 
spreading effect which occurs for smaller spatial dimen- 
sions: this spreading causes a reduction in the exchange 
splitting. At larger sizes, the exchange splitting is larger 
for the cuboidal shape due to the effect of the smaller 
height dimension of the cuboid, which pushes electron 
and hole together and thus increases the exchange split- 
ting. At larger dot sizes, the interaction scales like 



V. INTER-DOT INTERACTIONS 




FIG. 9: The intra-dot exchange Coulomb interaction strength 
integrand, plotted in the — Kz plane of reciprocal space, 
and around tire K = point. Note that it has value zero at 
K = and is asymmetric in x and z due to the choice of the 
Pz Kronig-Penney state for the holes. We have calculated the 
integrand for dot I and used the cubic geometry (a = /ii/2 = 
10 nm) and assumed that Ve = Vh = 500 meV. 




Dot Size (nm) 



In this section we shall discuss the m^er-dot coupling 
terms which are due to the Coulomb operator introduced 
in previous sections. These terms are crucial to the op- 
eration of a quantum device, since they may allow qubit- 
qubit interactions to take place, which is an essential 
requirement for two (or more) qubit gates to be con- 



FIG. 10: The intra-dot exchange Coulomb interaction 
strength, as a function of dot size and confinement po- 
tential. The solid lines are for a cubic geometry (a — hi/2), 
and the dotted Unes are for a cuboid (a = 5/ii). 
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structed. There are two important types of interaction 
which may occur. The first type is called the Forster 
interaction, and is described by an off-diagonal matrix 
element (in the computational basis) between two single 
exciton wavefunctions of the type introduced in Eqs. [S] 
andini but where the two excitons are located on different 
quantum dots (this interaction is called Vp in Eq.^. The 
second interaction which is important for this scheme is 
the direct self Coulomb interaction in a biexciton (double 
substitutional Slater determinant) wavefunction, where 
one exciton is located on each dot. This is called Vxx 
in Eq. ^ ^^^d amounts to the Coulomb binding energy 
between two excitons located on adjacent dots. We next 
quantify both of these interaction terms and discuss their 
properties within the context of the quantum computing 
implementation described in Sec. [n] . 



A. Off diagonal coupling: Forster interaction 




0.1 ^1.1 



The Forster or Coulomb exchange interaction can in- 
duce the transfer of an exciton from one quantum dot 
to the other. This is a non-radiative energy transfer 
whereby an exciton is destroyed on one dot and recreated 
on the other; it is an electrostatic interaction which pro- 
ceeds via a short lived virtual photon. Forster's original 
theorjiS^ showed that the interaction is dipole-dipole to 
lowest order; this theory was subsequently elaborated by 
Dexter>SSi who derived higher order and exchange terms 
in studies of the sensitized luminescence of solids. Here 
we extend this theory to the case of the many-body ex- 
citon states of quantum dots. The off-diagonal nature 
of the interaction causes the eigenstates of the Hamilto- 
nian, Eq.^ to become linear combinations of the compu- 
tational basis states |10) and |01). As described in Sec- 
tionm the degree of this mixing is crucial in determining 
how to generate quantum entanglement in the quantum 
dot molecule. The Forster coupling can be expressed as 
the matrix element of the direct Coulomb operator be- 
tween excitons located on each of the two dots: 



C / / iJs{VsWs(Vs) X 



(17) 



e^(R- 



rt)|R- 



■ ■(/'t(rt)V't(i"t)rfrsdrt 



This equation is equivalent to Eq. |3 but we have ex- 
plicitly included the interdot vector R in the Coulomb 
operator; we assume in this case that the two variables 
Ts and rt are defined from the centres of dot I and dot 
II respectively. We may evaluate Vp in exactly the same 
way as we evaluated the intra-dot couplings, so long as 
the new positions of the wavefunctions are included in 
the calculation. An example of the integrand appearing 
in Eq.Elis shown in Fig. ^2 It is interesting to compare 
this figure to Figs. 13 and |2| in the plot of Fig. [TTl there 
is an extra modulation due to the extra factor associated 
with the interdot separation, and this added modulation 
means the integral takes longer to evaluate numerically. 



FIG. 11: The plot of the integrand in K space which leads 
to the Forster strength, as a function of and {Ky = 
0). The extra modulation on the function (as compared with 
Figs. [7| and 1^ is caused by the interdot separation. The plot 
is for cubic dots with a = 6 = 10 nm and R — 20 nm, and we 
have taken Vs = Vh ~ 500 meV. 



The results are displayed in Figs. El and ^1 where the 
Forster strength is displayed as a function of dot sep- 
aration, shape and confining potential. The data are 
displayed on a log scale, and it can be seen that, for 
the cubic shape, they closely follow a 1/i?'^ law for all 
the separations considered. This form is expected for a 
dipole-dipole interaction, and we shall now discuss how 
a power series expansion and subsequent approximation 
leads to this type of interaction in this case. In so doing, 
we shall also explain why the interaction is modified as 
the size and shape of the dots are changed. 

By making the assumption that R is much larger than 
Ts and rt, we may Taylor expand the Coulomb operator. 
This procedure yields, to lowest non-zero order: 



Vf 



C 

7r3 



;ri) 



(r„)-— ((ri).R)((r„) 



R) , (18) 



where it has been assumed that the dielectric constant 
is independent of R + rg — rt , and as throughout the pa- 
per is assumed to take the constant value of = 10. 
The matrix element of the position operator between an 
electron and a hole state on dot I or II is 



(19) 



(i-i/ii) = / V^J/n(r)ri/;i/n(r)dr. 



Equation 1181 is therefore equivalent to the interaction of 
two point dipoles, one situated on each dot. We can 
proceed further by again employing the envelope func- 
tion approximation for electrons and holes (Eq. |3Jl and 
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FIG. 12: Dependence of the Forster interaction strength on 
the interdot separation 7?. The dots have equal sizes and re- 
sults are shown for two shapes: (i) cubic, with a = h/2 = 2 nm 
(upper two curves) and (ii) cuboidal, with a = 5h = 2 nm 
(lower two curves). The circles and squares represent the 
predictions of a full numerical simulation for well depths of 
500 meV and 2000 meV respectively (the electron and hole 
wells are assumed to be of the same depth). The dotted lines 
represent the predictions of the dipole-dipole model for the 
cubic shaped dots and the dashed curve represents the max- 
imum coupling predicted for the dipole-dipole model, which 
corresponds to Oi — On — 1.0. The dotted lines can be ob- 
tained by multiplying this maximum by the relevant values of 
Oi which can be obtained from Fig. 1141 



by rewriting Eq. ^| as 

(ri/„) - ^ V,,u I 0i/„(r - Ti) U'{i 

rrri 1 cell 



{Ti} 



(r-Ti)0i/„(r-Ti)[/(r)dr, 



(20) 



where Tj represents the set of lattice vectors. We have 
made use of the periodicity of the Bloch part of the wave- 
functions and assumed that this part of the wavefunction 
is the same for both dots. By making the assumption that 
the envelope function is slowly varying on the lengthscale 
of the atomic lattice and by using the orthogonality of the 
electron and hole Bloch functions wc find that: 



C 

~R3 



OiOii |(ra)| 



i?2 



((l-a) ■ R) 



(21) 



where the term (ra) represents the atomic position oper- 
ator expectation value 



(ra) 



Ue{v)vUh{v) dr, 



(22) 



FIG. 13: Dependence of the Forster interaction strength on 
the shape of the dots. The dots are assumed to be identical, 
but have a series of cuboidal shapes with different aspect ra- 
tios. The well depth is 500 meV (the electron and hole wells 
are assumed to be of the same depth) . The solid curve repre- 
sents the maximum coupling predicted for the dipole-dipole 
model, which corresponds to Oi — On — 1.0. 




Ve = Vh = 200 meV 

Ve = Vh = 500 meV 

— Ve = Vh = 1000 meV 

— - Ve = Vh = 2000 meV 
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FIG. 14: The overlap integral, Oi (Eq.ESJ, 

as a function of 

dot size and confinement potential. This graph can be used 
in conjunction with Fig. 1121 to obtain values of the Forster 
strength for a range of dot sizes and confinement potentials. 
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which is the same for both dots and 

Oi - / <?ii(r)<^LWdr, (23) 

space 

is the overlap of electron and hole envelope functions on 
the appropriate dot i. Eg. 1211 shows how the effects of 
the quantum dot size and shape (which determine the 
overlap integrals) may be separated from the effects of 
the material composition of the dot (which determine 
the atomic dipolc operator). 

It is now possible to obtain the strength of the inter- 
action by assuming the specific forms for the envelope 
and atomic functions which we discussed earlier, in Sec- 
tion lllll For the Kronig-Penney model, with a well width 
of 2x, the atomic position operator expectation value is 
given by 

(ra) = 32a;/9^2_ (34) 

The overlap integrals are easily calculated for the enve- 
lope functions described earlier, and are displayed as a 
function of dot size and confinement potential in Fig. E| 
We show the overlap integral for the usual two dot 
shapes: cubic (a = h/2) and flat cuboidal (a = 5/i), 
where the latter is more typical of dots grown by the 
Stranski-Krastanow method. The overlap is enhanced 
when there is a larger confinement potential and for 
larger dots, since in these cases the shape of the wave- 
function is less sensitive to the effective mass difference 
of the electrons and holes. We may use Fig. ^1 together 
with the atomic dipole value, to calculate the Forster 
strength for a range of dot sizes and confining potentials. 
Owing to its dependence on the atomic dipole operator 
iy-p oc x"^), we plot Vp/x'^ as a function of R in Fig. ^1 
Two example curves, for two equally sized cubic dots are 
shown for equal electron and hole potentials of 500 mcV 
and 2000 meV in Fig. ^1 together with the earlier full 
calculation in which the dipole-dipole approximation was 
not made. The full calculation was carried out for (i) cu- 
bic dots with a = h/2 = 2 nm, and (ii) cuboidal dots 
with a ^ 5h = 2 nm. For both shapes it is clear that the 
influence of dot shape and size is much more important in 
determining the size of the interaction than the influence 
of the size of the confinement potentials. Furthermore, 
the dipole-dipole approximation is very good in the case 
of a cubic dot, even at interdot separations which are 
relatively small when compared to the dot sizes. For 
the cuboidal dot, the dipole-dipole approximation fails 
at smaller separations. We investigated this effect more 
thoroughly by repeating the calculations for cuboids of 
different aspect ratio (see Fig. I13|l . The next order term 
(dipole-quadrupole) is zero in all cases when the dots 
have equal size; presumably the reason for the accuracy 
of the dipole-dipole approximation in the case of cubic 
dots is that the dipole-dipole terms dominate the higher 
order terms even at smaller dot separations. 

The simple Kronig-Penney model shows how the size 
of the Forster transfer depends upon the physical size 



of the atomic part of the wavefunction. However, (r'^) 
is a widely measured quantity since it determines the 
strength of dipole allowed transitions in optical spectra. 
In CdSe QDs it can be in the range of 0.9 to 5.2 eA,^^ in 
atomic systems it can also be several eAJ^ and in BSs and 
MSs has recently been observed to be about 1.7 eAiS 

As an illustration of the use of these curves, let us as- 
sume that we have a dot system in which, as before, ii = 
5 nm, a = 10 nm, 6 = 8 nm, and /ii = /i2 = 2 nm. Fur- 
thermore, let us take the measured dipole value for CdSe 
dots of 0.9 to 5.2 eA.'^^ In this case, the Forster strength 
is between 0.013 and 0.45 meV, which if Aq = would 
correspond to an on resonance energy transfer time of 
between 318 and 9.2 ps in dots with Vh = Ve = 500 meV. 
This is short enough to be useful for quantum computing 
purposes: decoherence times as long as a few na^ have 
been observed in QDs. In MSs or BSs, the interacting 
units can be as close together as 1 nm; using this and 
taking a typical molecular or biomolecular dipole value 
of about 1.7 eAfSL^ we obtain an interaction strength 
of 8.3 meV (or a transfer time of ~497 fs). Furthermore, 
Vf must certainly be controlled if the alternative scheme 
using Vxx is to be implemented (and therefore cannot 
be neglected as in Ref. Is^) . We note that Vp is not 
particularly sensitive to differences in dot size, though 
the differences in the diagonal (self energy) parts of the 
Hamiltonian which are caused by having dots of unequal 
size are very significant. We shall discuss this further in 
Section 

In the Section FV Bl we shall discuss how the biexciton 
binding energy term depends upon applied electric field, 
and how such a field may be crucial to the operation of 
a potential quantum logic device. We now discuss how 
the Forster term would vary when such a field is applied. 
Since an electric field would move the electron and hole 
away from one another, the overlap integrals, Eq. 
would be reduced by such a field. We show this specif- 
ically by simulating the effect of applying a field on the 
overlap integral Oi in Figs. 1151 andll6l In Fig. ^1 we see 
that Oi is significantly suppressed for fields of a few lO's 
of kV/cm, and that the suppression is easier to achieve 
in dots which have a larger dimension in the direction 
of the applied field. The reason for this is that in larger 
dots the electron and hole are more easily separated since 
there is more distance between the two dot-barrier inter- 
faces. In zero field, dots with a smaller dimension have 
smaller overlap integrals, for the reasons associated with 
the balance of kinetic and potential energy discussed ear- 
lier. Hence, we see that the curves for different dot sizes 
cross each other in an applied field. In Fig. ^| we plot 
the dependence of the Oi on the depth of the confinement 
potentials, for a = 10 nm for both the cubic and cuboidal 
geometries. We see here that the effect of varying con- 
finement potential is much smaller than varying dot size. 
The small difference that is evident, that of a slightly 
easier suppression for deeper confinement potentials, is 
presumably due to the fact that the wavefunctions in a 
shallower wells tend to be more spread out (since the 
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FIG. 15: The overlap integral, Oi (Ea. I23|l. as a function of 
electric field strength E for a range of dot sizes. The up- 
per figure shows the dependence when the dots take a cubic 
shape (a — hi/ 2); the lower figure shows the dependence for 
a cuboidal shape (a = 5/ii). Note that the overlap integral, 
and so also the Forster interaction, is suppressed at large field 
as the electron and hole are forced apart. 



potential energy cost in doing so is smaller), and so the 
overlap between an electron and hole at opposite sides of 
the dot is slightly enhanced. Again the curves cross one 
another since, as we discussed earlier, in zero field the de- 
localization of the states for shallower potentials means 
that the shape of each particle's wavefunction depends 
more strongly on its effective mass. 

The fact that the Forster coupling may be suppressed 
by an external field could be very useful: if an entangled 
state is produced by using this coupling, it may be main- 
tained by switching off Vp- If this could be done in a suf- 
ficiently short time (i.e. on the timescale of the evolution 
of the quantum device under the Forster coupling Hamil- 
tonian, but much less than typical decoherence times) , it 
may be possible to fabricate a two qubit gate using this 




E (kV/cm) 



FIG. 16: The overlap integral, Oi CEa. I23|l. as a function of 
electric field strength E for a range of confinement potentials, 
for a dot size of a = 10 nm. The upper figure shows the 
dependence when the dots take a cubic shape (a = h\/2); 
the lower figure shows the dependence for a cuboidal shape 
(a = 5/ii). 



effect. It may also be possible to achieve this switching 
in an alternative way, by leaving the Forster coupling at 
a constant value but by tuning the single exciton level 
spacing Aq through the electric field which induces the 
Stark shiftfi we shall return to this in Section IVTl ft is 
also interesting that a negligible Forster coupling is es- 
sential for the energy selective dot device discussed in 
Section ^ such a negligible coupling may be achieved 
through using an external field. There are disadvantages 
doing this in this case, however. A smaller electron-hole 
overlap will also reduce the coupling to the light field it- 
self, which we need to be strong enough to be able to 
perform conditional gates in a time short enough when 
compared with typical decoherence times. Hence, as so 
often in quantum computing implementation, a compro- 
mise must be struck between these two requirements. 
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B. Diagonal coupling 

We now calculate the direct Coulomb interaction be- 
tween two excitons, where one exciton is located on each 
dot. This interaction leads to the energy selectivity of 
the gate and is the responsible for the Vxx term of Eq. ^ 

Consider the following double substitutional Slater de- 
terminant, which we write using the hole prescription 
described in Section llll Bl It represents a combination of 
two ground conduction electron states and two ground 
hole states which correspond to one exciton on each dot: 



XX 



A 



i^l{r,)^Ur^)^l'{r, - RX(r4 - R) , (25) 



where A indicates that the wavefunction has overall an- 
tisymmetry, this being achieved by adding terms with 
labels swapped around in a Slater determinant form. R 
is the vector connecting the two dot centres, ri and ra 
represent the position vectors of electrons relative to the 
centres of dot I and dot II respectively and r2 and are 
the equivalent vectors for holes. The associated Coulomb 
operator Vxx is given by 



XX 



1 



IR + n-ral |R + ri-r4 
1 1 



|R + r2-r3| |R + r2-r4| 



(26) 



Expanding this expression in a Taylor series about R 
gives, to lowest non-zero order: 

Vxx = ^ jpi • Pii - ]|(Pi ■ R)(Pn • R)|, (27) 

where pi — e(ri — r2) is the overall dipole moment on dot 
I and pii — e(r3 — r4) is the overall dipole moment on dot 
II. To evaluate the matrix element {^xxl^xxl^'xx}, Pi 
and Pii in Ea. l27l are replaced by their expectation values 
for the wavefunction, Eq.lJHl This procedure gives rise to 
a direct term and exchange terms. The exchange terms 
arise from the parts of the wavefunction (Eq. I25|) which 
do not appear explicitly within the bracket but which 
have their labels swapped around and they are zero in 
the absence of wavefunction overlap between dots. The 
direct term is obtained through the use of the envelope 
function approximation, Eq. |3| which leads to the follow- 
ing equation for the expectation value (ri) 



(ri)= / ^i*(ri)ri0i(ri)dri 

J space 



(28) 



where the orthogonality of the Bloch functions for differ- 
ent bands and the slow variation approximation for the 
envelope functions have again been used. Similar expres- 
sions hold for the other position expectation values. 

For a cubic dot, where the electron and hole wavefunc- 
tions have a definite parity about the dot centre, Eq. [2H1 
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FIG. 17: (a) Exciton-exciton binding energy and (b) induced 
dipole moment as a function of the dot size, shape and applied 
electric field. We have assumed that 14 = = 500 meV. 



implies that the exciton-exciton coupling is zero. How- 
ever, this is not the case when this symmetry is broken. 
For instance, for pyramidal shapes the electron may lo- 
calize in one region of the dot and the hole in another 
regionii^ Alternatively, an electric field would induce a 
polarization on the dot; this field may be externally ap- 
plied or arise from intrinsic piezoelectric effects42iS 

We have simulated the effect of applying an electric 
field in our cuboidal model by including a linear poten- 
tial in the single particle Schrddinger Eqs. 01 and the 
results are displayed in Fig. El In the lower part of the 
figure we display the size of the exciton dipole moment 
Pi on one of the dots as a function of the dot size and of 
the applied electric field strength E. We do this for our 
usual two geometries: cubic (a = hi/ 2) and flat cuboid 
(a = bhi) (see Fig. 0J). In both cases the field is apphed 
along one of the axes of the square base. The interac- 
tion strength Vxx is then obtained by using the size of 
Pi for each dot and substituting into Eq. |^ Thus, the 
upper part of Fig. 1171 shows the strength Vxx, normal- 
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ized by R^, for two dots of equal size and calculated for 
both of the dot geometries just described. (The interac- 
tion between two unequally sized dots can similarly be 
obtained by the use of Fig. [TTJb) and Eq. [23) At very 
small applied field, the induced dipolc is linearly pro- 
portional to the field, and hence the interaction strength 
takes a quadratic dependence on field. At larger applied 
fields, the induced dipole begins to saturate as the elec- 
tron and hole approach the edges of confining potential 
of the quantum dots; this limits the useful interaction 
strength which may be obtained from a given pair of 
dots. It is interesting that the interaction strength is 
much more dependent on a than it is on hi; this is be- 
cause it is in the basal plane that the field is applied, 
and so it is in this direction that the dipole moment is 
induced. The relative insensitivity of Vxx to hi turns 
out to be very useful: it means that SK dots, which can 
be stacked closely on top of one another but which have 
a relatively large base size, can be made to interact very 
strongly. As an illustration, consider the following typi- 
cal SK dot parameters: i? = 5 nm, a = 8 nm, 5 = 10 nm, 
and hi — h2 = 2 nm (we have assumed the usual experi- 
mental situation in which the upper dot of the stack has a 
slightly larger size). In an x directed field of 100 kV/cm, 
these parameters give Vxx ~ 120 meV, which would re- 
sult in a lower time limit for the gate operation of around 
10 fs. This is relatively short; decoherence times on the 
order of nanoseconds have been observed for uncoupled 
dotsiStt Finally, we note that we have only calculated Vxx 
to first order; in some cases higher order terms may be 
important. 




VI. FURTHER DISCUSSION 



The model outlined above can also be used to calcu- 
late Aq, the difference in exciton creation energy for two 
different sized dots in the absence of interactions. This 
is done by simply calculating the single particle electron 
and hole energy and taking into account the Coulomb 
binding energy between them. We also assume that the 
electron and hole are in the spin singlet state; spin is not 
important for our present proposal and it is always possi- 
ble to choose the spin singlet state by using light with the 
appropriate polarization. Hence we can effectively ignore 
the triplet states in considering the scheme described in 
Section ini 

The absolute value of Aq is displayed as a function of 
the ratio of the dot side lengths for each of the usual 
geometries in Fig. ^1 We have displayed Aq, both in 
the absence of the (intra-dot) Coulomb terms, and when 
these terms are included to first order by using the cal- 
culations of Section IIVI Aq is zero when the dots are 
of equal size, and then increases as the difference in size 
becomes greater. The Coulomb terms serve to reduce 
the size of Aq because they are larger for the smaller 
dot of the pair (which also has the larger single particle 
energies). 



FIG. 18: Energy splitting Ao = uji—uj2 of the singlet qubit ex- 
citon states l^oi), and I'I'io) for (a) cubic dots (a = hi/2, h = 
h2/2) and (b) cuboidal dots (a = 5hi,b = 6^2) in the ab- 
sence of the Forster interaction as a function of the dot size 
ratio a/b. The splitting is independent of interdot distance. 
The solid lines represent the splitting in the presence of the 
Coulomb and exchange splitting terms and each adjacent dot- 
ted line represents the splitting without Coulomb interactions 
for each dot size. 



We now use the analysis of Section [H] and the calcula- 
tions of subsequent sections to obtain the size of the ci 
component of the l^'io) and |4'oi) states. If we assume 
that the Forster strength is small in comparison with Aq, 
then we have that ci w Vp/Aq (see inset of Fig. ^c)) and 
by substituting Eq. |31 and Eq. 1211 we obtain 



R^ci 



37.1 



x^OiOu Ao 



(29) 



where Aq is measured in meV, and x and R are in nm. 
This quantity is displayed as a function of dot size ratio 
for the usual geometries in Fig. E| It can be seen there 
that a range of ci values can be obtained by choosing dots 
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FIG. 19: The size of the component of the wavefunction ci 
for (a) cubic dots (a = b = /i2/2) and (b) cuboidal dots 

(o = 5hi,b = 5/12) as a function of the dot size ratio a/b. ci 
has been scaled by its dependence on the interdot distance, 
R, typical atomic spacing, x, and overlap integrals, Oi. 



with appropriate values of a;, R and a/b. For example, 
cubic dots with large x (> 1 nm say), small R {< 3 nm 
say) and a/& ~ 1 give a larger ci, and it is then more 
appropriate to use the Forster interaction itself to create 
entangled states. On the other hand, dots with smaller x, 
larger i? or a large mismatch in dot size would be more 
suited to the scheme which uses the T^x for QC and 
entanglement generation. A scheme similar to the latter 
one was discussed by Biolatti et al. in Ref . 53, though the 
ofF-diagonal coupling was not considered there at all; we 
now see how important it is to consider the effect of this 
interaction. The fidelity^^ of a typical Vxx entangling 
gate operation (e.g. |11) |10)) is equal to 1 — — and 
so one must be careful when using the biexciton scheme 
to use the available parameter space and make sure that 
the Forster transfer is suppressed to the desired accuracy. 
There are other sources of decoherence in this case (e.g. 



the interaction with optical and acoustic phononsSiSLSS) 
which will reduce the value of the fidelity to below 1 — cf . 
To minimize the effects of such decoherence channels, it 
is important to maximize Vxx, since this leads to an im- 
proved transition discrimination and so to a faster gating 
time. This can be done by applying an electric field and 
choosing an appropriate dot shape, size and separation 
(as described earlier). It is then necessary to minimize 
the basis state mixing for the chosen parameters by se- 
lecting a suitable dot size ratio and material composition. 
It was seen earlier that the value of Vxx could be as high 
as several tens of meV. If we assume a conservative fig- 
ure of 10 meV, we find that the uncertainty principle 
implies that a CNOT gate could be performed in a time 
of around 100 fs. This is relatively short; decoherence 
times on the order of nanoseconds have been observed 
recently for uncoupled dots.— Hence we conclude that 
this scheme looks rather promising as a solid state im- 
plementation of quantum computation. 

It is also possible to use the Stark effeclii to tune two 
non-resonant levels into resonance, thereby allowing for 
the kind of switching of the Forster interaction which is 
required if it is to be used for quantum logic. This can 
be done so long as the two dots arc made such that they 
have different polarizibilities (which could be achieved 
by using either different sized dots or dots made of dif- 
ferent materials). So long as this difference is such that 
the levels are brought closer together by applying a field, 
and that shifts as large as Aq can be achieved, switch- 
able resonant transfer is possible. It might be difficult to 
achieve the switching in a time which is short enough for 
a quantum gate to be performed; however, in this case 
the optical (AC) Stark shift could be employed by using 
ultrafast lasers iSS 

Single shot qubit state measurement in QDs could 
be performed by using resonant fluorescent shelving 
techniquesiS The QD state measurement can also be 
achieved by means of projecting onto the computa- 
tional basis and measuring the final register state by ex- 
ploiting ultrafast near-field optical spectroscopy and mi- 
croscopyiS2i2i these allow one to address, to excite and 
to probe the QD excitonic states with spectral and spa- 
tial selectivity. In addition, the qubit register density 
matrix can be reconstructed by measuring the QD pho- 
ton correlations via standard quantum state tomography 
techniques. In particular, we believe that the activity 
of the (Forster) resonant energy transfer processes dis- 
cussed in this paper can be accomplished in our coupled 
dot molecule by measuring the intensity correlation func- 
tion (usually denoted g*-^'') in a Hanbury-Brown/Twiss 



type experiment 



6.8.86 



Such an experiment can reveal sig- 



natures of purely non-classical photon correlations aris- 
ing from the QD molecule emission (i.e., photon anti- 
bunching or bunching behaviour). This idea has already 
been experimentally explored for the case of pairs of dye 
molecules by Berglund et al.S^ Scalability of the scheme 
given here could also be possible by adopting a globally 
addressed qubit strategySl on a stack of self-organized 
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We finalize this section with a discussion on how the 
Vp coupling can be usefully manipulated in biomolecular 
nanostructures. Light-harvesting antenna complexesi^ 
or arrays of strongly interacting individual molecules^ 
could provide an appropriate system in which the Forster 
interaction could be used for QIP tasks. They are gener- 
ally very uniform structures, and we may compare them 
to QDs by setting a/6 ~ 1, or Vf/Aq ;» 1. Then the one- 
exciton eigenstates of a two qubit system with a Forster 
coupling naturally allows the generation of the states 
^(|01) ± |10)), which, apart from their applications to 
quantum protocols, can be particularly useful in the fight 
against decoherence. Spectroscopic, line-narrowing tech- 
niques (e.g., hole burning and site-selective fluorescence), 
infrared and Raman experimental studies reveal that the 
main decoherence mechanisms in the antenna complexes 
arise from energetic disorder, electron-phonon coupling, 
and temperature effectsi^ In this scenario, the excita- 
tions couple to an environment that typically possesses 
a much larger coherence length than the biomolecular 
units (BChl's) spacing. For example, the BChl's in the 
antenna complex LH2, which we regard as a potential 
system for quantum logic, are spaced by as little as 1 nm, 
and hence so-called collective decoherence is expected to 
apply. In this case, provided that the logical qubit en- 
coding li).^ = \01)^^., IT), = |10)^.j. that uses two phys- 
ical (exciton) qubits can be realized in the BChl's sys- 
tem, arbitrary superpositions of logical qubits such as 
{a^ li). + A IT),)^"^, i = 1, . . . , Af, ft e C, are im- 
mune to dephasing noise (described by a tr^ operator—), 
and single qubit manipulations can be carried out on the 
timescale of the Forster coupling (which as we have seen 
can be as short as 497 fs). Two-qubit logic gates can 
also be implemented within a decoherence-free subspace 
by using the above encoding)2£ thus completing a uni- 
versal set of gates. Initialization of the system requires 
the pairing of the physical qubits to the logical 'ground' 
state |i)f^, and readout is to be accomplished by iden- 
tifying on which of the two structures the exciton re- 
sides. Furthermore, rings of BChl's appear side by side 
in naturally occuring antenna complexes and also display 
energy selectivity — smaller rings tend to have higher en- 
ergy transitionsi^ Thus, following a scheme as above, 
it may be possible to scale up such biological units in 
a natural way and construct a robust energy selective 
scheme for quantum computation.^^ We also note that 
arrays of strongly interacting individual molecules which 
are coupled via a near field dipole-dipole interaction are 
well suited for our quantum computing and entangling 
schemes, especially due to the existence of Vxx-type of 
energy shifts, which have been analyzed in this paper. 
A full discussion of these ideas will be presented else- 
where 



VII. SUMMARY 

We have shown that it is possible to use the two dif- 
ferent electrostatic coupling terms (Forster transfer and 
biexciton binding energy) between excitons in quantum 
dots to construct two qubit gates which, in addition to 
an appropriate control over single qubits, are enough for 
universal quantum computation. We have also discussed 
how to generate tailored exciton entangled states by us- 
ing these gates. We have furthermore modelled a pair of 
quantum dots in the simplest envelope function approx- 
imation, and have mapped out the areas of parameters 
space where one of the two interactions dominate. We 
have discussed in detail how to perform entangling oper- 
ations when one of the two interactions is dominant; the 
case when the two have similar magnitude leads to a rich 
spectrum of entangled states whose degree of entangle- 
ment can be quantified a posteriori in, for example, pho- 
ton correlation experiments. We have concentrated here 
on two geometries, namely that of a cube-shaped dot and 
that of a cuboid-shaped dot, and our calculations have 
been partially analytical and partially numerical. In fu- 
ture, we are hoping to obtain simpler analytical results 
for many of the quantities we have calculated, by using 
different dot geometries. We shall also extend our calcu- 
lations to include decoherence effects, and in future work 
we shall be particularly interested on how we might use 
the Forster interaction to create and use decoherence free 
subspaces. 
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APPENDIX A: SINGLE PARTICLE SOLUTIONS 

As we discussed in the text, the solution of the three 
dimensional finite well square box potential is obtained 
by using the solutions of the one dimensional finite square 
well and expanding the Schrodinger equation in these ba- 
sis states. In order to do this, both bound and unbound 
basis states had to be taken into consideration. 



1. Bound States 

The problem of finding the solutions of a finite square 
well is covered in most undergraduate text books (see 
e.g. Ref . WH) . and so we shall not go into too much detail 
here. If we assume that the finite well is centred around 
x = 0, we have different forms of solution depending on 
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FIG. 20: Schematic diagram of the potential used to generate 
the unbound basis states used in the calculations. 



the parity of the wavefunction. With reference to Fig. 1201 
and assuming that our wavefunctions have decayed once 
the infinite barriers are reached (see Section IIII A|l , we 
have 



Region 


Even Solutions 


Odd Solutions 


X < -a/2 
-a/2 < X < a/2 
X > a/2 


B cos{kx) 


D cos{kx) 



These equations were solved by using the numerical root 
finding algorithm provided with the NAG package. Once 
the energies are obtained, a and k follow from simple 
substitution. The normalization constants follow from 
the usual normalization procedures once all of the other 
parameters are known. 

2. Unbound States 

The problem for unbound states is somewhat less 
straightforward. We assume that the finite well (of width 
2a) is embedded within an infinite well (of width 2L), and 
that we may set the wavefunction outside the infinite well 
to zero (see Fig. 12011 . Then there are again three regions 
which have different forms of solution, which take the 
form: 



Region 


Even Solutions 


Odd Solutions 


-L/2 <x < -a/2 
-a/2 < X < a/2 
a/2<x < L/2 


-Asm{k'{x- L)) 

B cos{kx) 
Asui{k' {x - L)) 


C8m{k'{x - L)) 

D cos{kx) 
-C sm{k' {x - L)) 



TABLE II: Table of the forms of the unbound state solutions 
for the finite square well potential. 



The wavevector inside the dots is given by fc = 
{2m*E)i/h, and k' = {2m* {E - Vp))i/h. The solution 
is obtained by again invoking the continuity of the wave- 
function amplitude and the probability current at the 
boundaries, and the following transcendental equations 
are found for the even and odd solutions respectively: 



TABLE I: Table of the forms of the bound state solutions for 
the finite square well potential. 



E-Vp 
E 



/ a{2m*E)2 \ 
tan — ; — — X 



(A3) 



Where A, B, C and D are normalization constants. 
a — {2mp{Vp— E))^ / h and k = {2m*pE)2 /h. By ensuring 
the continuity of the amplitudes of the wavefunctions and 
the probability current at the boundaries, the following 
transcendental equations for the energy E are obtained: 

For even solutions: 



tan 



{2m;EY 



2h 



For odd solutions: 



Vp~E 
E 



(Al) 



E-Vp 



tan 



{a-L)i2m,*{E-Vp))2 



cot °(^-:^)Mx 



(A4) 



tan 



[a - L){2m*{E - Fp))2 

h 



cot 



{2m;E)2a 
2h 



Vp-E 
E 



Eas. lA3l and lA4l are solved for E, and then the wavevec- 
(A2) tors and normalization constants follow on as before. 



* Electronic address: Ibrendon.lovettQmaterials. oxford. ac.uki 

t On leave of absence from Centro Internacional de Fi'sica 
(GIF), A. A. 4948, Bogota, Colombia; Electronic address: 
j.reina-estupinanOphysics. oxford. ac.uk 

^ P. Harrison, Quantum Wells, Wires and Dots (Wiley, New 



York, 2001). 

^ D. Bimberg, M. Grundmann, and N. N. Ledentsov, Quan- 
tum Dot Hetero structures, (Wiley, Ghichester, 1999). 

^ M. Bayer, O. Stern, P. Hawrylak, S. Fafard, and A. Forchel, 
Nature 405, 923 (2000). 



19 



* H. Kamada, H. Gotoh, J. Tommyo, T. Takagahara, and 
H. Ando, Phys. Rev. Lett. 87, 246401 (2001); H. Htoon, 
T. Takagahara, D. Kulik, O. Baklcnov, A. L. Holmes 
and C. K. Shih, Phys. Rev. Lett. 88, 087401 (2001); 
T. H. Stievater, X. Q. Li, D. G. Steel, D. Gammon, 
D. S. Katzer, D. Park, C. Picrmarocchi, L. J. Sham, 
Phys. Rev. Lett 87 133603 (2001). 

^ P. Michlor, A. Imamoglu, M. D. Mason, P. J. Carson, G. F. 

Strouse, and S. K. Buratto, Nature 406, 968 (2000). 
^ D. Regelman, U. Mizrahi, D. Gorshoni, E. Ehrcnfrcund, 

W. V. Schoenfeld, and P. M. Petroff, Phys. Rev. Lett. 87, 

257401 (2001). 

P. Michler, C. Becher, W. V. Scoenfeld, P. M. Petroff, L. D. 
Zhang, E. Hu, and A. Imamoglu, Science 290, 2282 (2000). 

* C. Santori, M. Pelton, G. Solomon, Y. Dale, and Y. Yar- 
rnamoto, Phys. Rev. Lett. 86, 1502 (2001). 

J. H. Reina, D. Phil. Thesis, Oxford University (2002). 
J. H. Reina, L. Quiroga, and N. F. Johnson, Phys. Rev. A 
62, 012305 (2000); L. Quiroga and N. F. Johnson, 
Phys. Rev. Lett. 83, 2270 (1999). 
" A. M. Steane, Rep. Prog. Phys. 61, 117 (1998). 

M. A. Nielsen and I. L. Chuang, Quantum Computation 
and Quantum Information (CUP, Cambridge, 2000). 
P. W. Shor, in Proceedings of the 35th Annual Symposium 
on the Foundations of Computer Science (IEEE Press, New 
York, 1994). 

" R. P. Feynman, Int. J. Thcor. Phys. 21, 467 (1982). 

D. S. Abrams and S. Lloyd, Phys. Rev. Lett. 79, 2586 
(1997). 

^® C. Zalka, Proc. R. Soc. London A 454, 313 (1998). 

See e.g., the issue on Practical Realizations of Quantum 

Information Processing, Phil. Trans. R. Soc. Lond. A 361, 

No. 1808 (2003). 
^® J. I. Cirac and P. ZoUcr, Phys. Rev. Lett. 74, 4091 (1995). 
1^ J. I. Cirac and P. Zoller, Nature 404, 579 (2000). 
^° C. Monroe, D. M. Meekhof, B. E. King, W. M. Itano, and 

D. J. Wineland, Phys. Rev. Lett. 75, 4714 (1995). 

K. M0lmer and A. S0rensen, Phys. Rev. Lett. 82, 1835 

(1999). 

C. A. Sackctt, D. Kielpinski, B. E. King, C. Langer, 
V. Meyer, C. J. Myatt, M. Rowe, Q. A. Turchettc, W. M. 
Itano, D. J. Wineland, et al.. Nature 404, 256 (2000). 
T. Pellizzari, S. Gardiner, J. Cirac, and P. Zoller, Phys. 
Rev. Lett. 75, 3788 (1995). 

Q. A. Turchettc, C. J. Hood, W. Langc, H. Mabuclii, and 

H. J. Kimble, Phys. Rev. Lett. 75, 4710 (1995). 

J. I. Cirac, T. Pellizari, and P. Zoller, Science 273, 1207 

(1996) . 

^® J. I. Cirac, P. Zoller, H. J. Kimble, and H. Mabuchi, Phys. 
Rev. Lett. 78, 3221 (1997). 

A. Rauschenbcutel, G. Nogucs, S. Osnaghi, P. Bertet, 
M. Brune, J. M. Raimond, and S. Haroche, Phys. Rev. 
Lett. 83, 5166 (1999); Phys. Rev. A 64, 050301(R) (2001). 
2* N. A. Gershenfeld and I. L. Chuang, Science 275, 350 

(1997) . 

I. L. Chuang, L. M. K. Vandersypen, X. L. Zhou, D. W. 
Leung, and S. Lloyd, Nature 393, 143 (1998). 

^° D. G. Cory, A. F. Fahmy, and T. F. Havel, Proc. Natn. 
Acad. Sci. USA 94, 1634 (1997). 

J. A. Jones and M. Mosca, J. Chem. Phys. 109, 1648 

(1998) . 

32 B. E. Kane, Nature 393, 133 (1998). 

33 R. Vrijcn, E. Yablonovitch, K. Wang, H. W. Jiang, A. Ba- 
landin, V. Roychowdhury, T. Mor, and D. P. DiVincenzo, 



Phys. Rev. A 62, 012306 (2000). 
3* D. Mozyrsky, V. Privman, and M. L. Glasscr, Phys. Rev. 

Lett. 86, 5112 (2001). 
35 G. K. Brennen, C. M. Caves, P. S. Jessen, and I. H. 

Deutsch, Phys. Rev. Lett. 82, 1060 (1999). 
3® D. Jaksch, H.-J. Briegel, J. I. Cirac, C. W. Gardiner, and 

P. Zoller, Phys. Rev. Lett. 82, 1975 (1999). 
3'' M. Greiner, O. Mandel, T. Esshnger, T. W. Hansch, and 

I. Bloch, Nature 415, 39 (2002). 
3* D. V. Averin, Solid State Commun. 105, 659 (1998). 
3^5 Y. Makhlin, G. Schon, and A. Shnirman, Nature 398, 305 

(1999) ; Rev. Mod. Phys., in press (2001). 

*° Y. Nakamura, Y. A. Pashkin, and J. S. Tsai, Nature 398, 
786 (1999). 

'^^ C. H. van dor Wal, A. C. J. tcr Haar, F. K. Wilhclm, 
R. N. Schouten, C. J. P. M. Harmans, S. Lloyd, and J. E. 
Mooij, Science 290, 773 (2000); J. E. Mooij, T. P. Orlando, 
L. Levitov, L. Tian, C. H. van der Wal and S. Lloyd, Sci- 
ence 285, 1036 (1999). 

^'^ A. Barenco, D. Deutsch, A. Ekert, and R. Jozsa, Phys. 
Rev. Lett. 74, 4083 (1995). 

"3 G. Burkard, D. Loss, and D. P. DiVincenzo, Phys. Rev. B 
59, 2070 (1999). 

J. H. Reina, L. Quiroga, and N. F. Johnson, Phys. Rev. B 
62, 2267(R) (2000). 

E. Biolatti, R. C. lotti, P. Zanardi, and F. Rossi, Phys. 
Rev. Lett. 85, 5647 (2000). 

B. W. Lovett, J. H. Reina, A. Nazir, B. Kothari, and 
G. A. D. Briggs, Phys. Lett. A 315, 136 (2003). 
S. de Rinaldis, I. D'Amico, E. Biolatti, R. Rinaldi, R. Cin- 
golani, and F. Rossi, Phys. Rev. B 65, 081309 (2002). 
X. D. Hu and S. Das Sarma, Phys. Rev. A 61, 062301 

(2000) . 

G. Chen, N. H. Bonadco, D. G. Steel, D. Ganmion, D. S. 

Katzer, D. Park, and L. J. Sham, Science 289, 1906 (2000). 
5° X. Q. Li, Y. W. Wu, D. Steel, D. Gammon, T. H. Stievater, 

D. S. Katzer, D. Park, C. Piermarocchi, and L. J. Sham, 

Science 301, 809 (2003). 
'^^ F. Troiani, U. Hohenester, and E. Molinari, Phys. Rev. B 

62, R2263 (2000). 

52 E. Pazy, E. Biolatti, T. Calarco, I. D'Amico, P. Zanardi, 

F. Rossi, and P. Zoller, Europhys. Lett. 62, 175 (2003). 

53 E. Biolatti, I. D'Amico, P. Zanardi, and F. Rossi, 
Phys. Rev. B 65, 075306 (2002). 

•"^ T. Forstcr, Disc. Farad. Soc. 27, 7 (1959). 

55 D. L. Dexter, J. Chem. Phys. 21, 836 (1953). 

5® X. Hu, T. Ritz, A. Damjanovic, F. Autenrieth, and 

K. Schulten, Quarterly Rev. of Biophysics 35, 1 (2002). 
5^" S. A. Crooker, J. A. Holhngsworth, S. Trctiak, and V. I. 

Klimov, Phys. Rev. Lett. 89, 186802 (2002). 
5* C. Hettich, C. Schmitt, J. Zitzmann, S. Kuhn, I. Gerhardt, 

and V. Sandoghdar, Science 298, 385 (2002). 
5^ P. Borri, W. Langbein, S. Schneider, U. Woggon, R. L. 

Sellin, D. Ouyang, and D. Bimberg, Phys. Rev. B 66, 

081306(R) (2002). 
®° B. Schelpe, A. Kent, W. J. Munro, and T. P. Spiller, 

Phys. Rev. A 67, 052316 (2003). 
®^ J. H. Reina, L. Quiroga, and N. F. Johnson, Phys. Rev. A 

65, 032326 (2002). 
®2 A. Nazir, B. W. Lovett, J. H. Reina, and G. A. D. Briggs, 

Preprint xxx. lanl.gov/quant-ph/0309099. 
®3 D. J. Eaglesham and M. Cerullo, Phys. Rev. Lett. 64, 1943 

(1990). 

Q. Xie, A. Madhukar, P. Chen, and N. P. Kobayashi, 



20 



Phys. Rev. Lett. 75, 2542 (1995). 

J. Tcrsoff, C. Teichert, and M. G. Lagally, Phys. Rev. Lett. 

76, 1675 (1996). 
'^'5 G. S. Solomon, J. A. Trezza, A. F. Marshall, and J. S. 

Harris, Phys. Rev. Lett. 76, 952 (1996). 
^"^ A. Franceschetti and A. Zunger, Phys. Rev. Lett 78, 915 

(1997). 

^® A. Franceschetti, H. Fu, L. W. Wang, and A. Zunger, 

Phys. Rev. B 60, 1819 (1999); ibid. 58, 13367(R) (1998). 
'^^ A. J. Williamson, L. W. Wang, and A. Zunger, Phys. Rev. 

B 62, 12963 (2000). 
™ H. T. Johnson, L. B. Freund, C. D. Akyuz, and A. Za- 

slavsky, J. Appl. Phys. 84, 3714 (1998). 

M. A. Cusack, P. R. Briddon, and M. Jaros, Phys. Rev. B 

54, R2300 (1996). 
''^ C. Pryor, M.-E. Pistol, and L. Samuelson, Phys. Rev. B. 

56, 10404 (1997); C. Pryor, Phys. Rev. B. 60, 2869 (1999). 
'^^ M. CaUfano and P. Harrison, J. Appl. Phys. 88, 5870 

(2000); ibid. 86, 5054 (1999). 
'''' S. Gangopadhyay and B. R. Nag, J. Appl. Phys. 81, 7885 

(1997). 

G. Bastard, Phys. Rev. B 24, 5693 (1981). 

^® A. Szabo and N. S. Ostlund, Modem Quantum Chemistry 
(McGraw-Hill, New York, 1989). 

"^""^ C. Cohen- Tanoudji, B. Diu, and F. Laloe, Quantum Me- 
chanics (Wiley, London, 1977). 

'■^ T. Forster, Ann. Phys. (Leipzig) 2, 55 (1948). 



™ http:/ /physics. nist.gov/Pubs/AtSpoc/nodel7.html. 

D. Birkodal, K. Loosson, and J. M. Hvam, Phys. Rev. Lett. 
87, 227401 (2001); P. Borri, W. Langboin, S. Schneider, 
U. Woggon, R. L. Sellin, D. Ouyang and D. Bimberg, 
Phys. Rev. Lett. 87, 157401 (2001); M. Bayer and A. 
Forchel, Phys. Rev. B. 65, 041308(R) (2002). 

®^ J. F. Poyatos, J. I. Cirac, and P. Zoller, Phys. Rev. Lett. 
78, 390 (1997). 

®^ C. Cohen- Tanoudji, J. Dupont-Roc, and G. Grynberg, 
Atom-Photon Interactions: Basic Processes and Applica- 
tions (Wiley, New York, 1992). 

®^ C. Roos, T. Zeiger, H. Rohde, H. C. Nagerl, J. Eschner, 
D. Leibfricd, F. Schmidt-Kaler, and R. Blatt, Phys. Rev. 
Lett. 83, 4713 (1999). 

J. R. Guest, T. H. Stiovatcr, G. Chen, E. A. Tabak, B. G. 

Orr, D. G. Steel, D. Gammon, and D. S. Katzer, Science 

293, 2224 (2001). 
®^ A. G. White, D. F. V. James, W. J. Munro, and P. G. 

Kwiat, Phys. Rev. A 65, 012301 (2002). 
*® A. J. Berglund, A. C. Doherty, and H. Mabuchi, Phys. 

Rev. Lett. 89, 068101 (2001). 

S. C. Benjamin, Phys. Rev. A 61, 020301 (2000). 
®® D. A. Lidar, I. L. Chuang, and K. B. Whaley, Phys. Rev. 

Lett. 81, 2594 (1998). 
®^ J. H. Reina, B. W. Lovett, A. Nazir, and G. A. D. Briggs, 

in preparation. 



